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Abstract 


A theory is developed for predicting wing-rock characteristics. 

From available data, it can be concluded that wing rock is triggered by 
flow asymmetries , developed by negative or weakly positive roll damping, 
and sustained by nonlinear aerodynamic roll damping. A new nonlinear 
aerodynamic model that includes all essential aerodynamic nonlinearities 
is developed. The Beecham-Titchener method is applied to obtain 
approximate analytic solutions for the amplitude and frequency of the 
limit cycle based on the three-degree-of-freedom equations of motion. 

An iterative scheme is developed to calculate the average aerodynamic 
derivatives and dynamic characteristics at limit-cycle conditions. Good 
agreement between theoretical and experimental results is obtained. 
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1 . 


Introduction 


Modem combat airplanes and tactical missiles have been designed to 
fly in the regime of high angles of attack to achieve higher maneuver- 
ability and a larger flight envelope. Flying at high angles of attack 
will often produce lateral-directional instabilities such as wing rock, 
wing drop, nose wander, and nose slice. The main topic here is wing 
rock. 

In this study, a summary based on extensive data is presented to 
explain the basic mechanisms of wing rock. Nonlinear aerodynamic models 
for one- and three-degree-of-freedom (DOF) motions are developed to 
investigate the main aerodynamic nonlinearities causing wing rock. The 
Beecham-Titchener method (ref. 1) is applied and extended to determine 
the approximate analytic solutions for the one- and three-DOF equations 
of motion. Because the interaction between aerodynamics and dynamics of 
wing rock is strongly nonlinear, an iterative scheme is developed to 
obtain the average aerodynamic derivatives and limit-cycle character- 
istics. 


1 



2 


Wing Rock Phenomenon 


Wing rock is an uncommanded roll-yaw oscillation dominated by roll 
motion oscillating with a constant amplitude (fig. 1). Most flight 
records show that the amplitude of the limit cycle in roll is at least 
one order of magnitude higher than that of the oscillation in yaw or 
pitch (fig. 2). Since wing rock is mainly a result of aerodynamic 
nonlinearities at high angles of attack, the characteristic motion must 
be described through a nonlinear mathematical aerodynamic model. 

However, the response of a lightly damped Dutch roll can be adequately 
determined by a linear mathematical aerodynamic model. Therefore, both 
mathematically and physically, wing rock is significantly different from 
a lightly damped Dutch roll. 

The severity of wing rock depends on the amplitude and period of 
the limit cycle. It may degrade weapon aiming accuracy (ref. 2), 
missile avoidance capability, and turning effectiveness. It may also 
cause safety problems during a landing approach or during a dogfight. 

To improve the air combat performance or handling qualities, the basic 
mechanisms of wing rock need to be understood. 

Wing rock strongly depends on the details of the configuration 
geometry. It is caused by the intricate flow pattern (ref. 3) at high 
angles of attack that exists around an aircraft, especially the wing, 
and is very sensitive to small changes in aircraft geometry. For 
example, an F-4 when installed with leading edge slats (ref. 4) showed 
its onset of wing rock deferred to higher angles of attack and the 
amplitude buildup became much milder and more gradual, in contrast to a 
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basic configuration. The shark nose and wing-root leading-edge 
extension (ref. 5) on an F-5 were very effective in suppressing wing 
rock. As compared with the Gnat aircraft carrying no fuel tanks 
(fig. 3) those carrying fuel tanks closely under its wings (ref. 6) 
exhibited wing-rock onset at much higher angles of attack. The early 
versions of the Harrier (ref. 7) equipped with two close fences on the 
wing suffered from unacceptable wing rock; but a modified Harrier 
equipped with two well-spaced fences flew with excellent roll steadi- 
ness. Recent wind-tunnel results (fig. 1) disclosed that a delta wing 
with 80-degree leading-edge sweep (refs. 8 and 9) would indulge in self- 
excited wing rock. On the other hand, another delta wing with 76-degree 
leading-edge sweep (ref. 9) did not exhibit such motion at all. 

Since the induced increment in angle of attack on a local wing 
panel varies periodically during wing rock, the flow pattern differs at 
different instants. For a moderately swept wing, unsteady turbulent 
boundary layers may alternately separate and reattach on the wing due to 
unsteady shock movement and boundary layer interaction. For a slender 
wing, leading-edge vortices may periodically vary their lateral and 
vertical positions (fig. 4). Furthermore, asymmetric forebody vortices 
(fig. 5) may interact with the wing in a nonsteady manner. As a result, 
the pressure field changes cyclically during wing rock. This in turn 
drives the aircraft to rock back and forth. Hence, the nonlinear inter- 
action between the unsteady flow field and the aircraft motion is 
extremely important for wing rock. 

Success of any attempt to control wing rock depends on the 
understanding of its dynamic and aerodynamic mechanisms. A unified 
theory is developed here to interpret the basic mechanisms of wing 
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rock. After examining a large amount of data, it is concluded that wing 
rock can be triggered by some asymmetric flow conditions, developed by 
the loss of roll damping, and sustained primarily by nonlinear roll 
damping. These are explained in the following. 

2.1. Triggering Mechanisms 

Wing rock may be initiated either during an asymmetric flight with 
a finite sideslip or during a zero-sideslip flight with some flow 
asymmetries over the aircraft. 

An airplane in asymmetric flight at high angles of attack, such as 

uncoordinated turns, asymmetric firing of weapons, and other maneuver 

situations, is susceptible to the initiation of wing rock. In addition, 

poor lateral-directional stability and ineffective control power may 

also trigger wing rock. For example, wing rock may be triggered due to 

(1) rapid roll-yaw control input for the Gnat aircraft (ref. 10) and the 

F-5E (ref. 11), (2) zero Dutch roll damping of the Gnat aircraft (ref. 

10) and the Handley Page 115 research aircraft (ref. 12), and (3) 

vanishing dynamic directional departure parameter (C ) of the F-4 

n 0 . 

(refs. 4 and 13) and the F-5 (refs. 5 and 14). 

Asymmetric flow occurring at zero-sideslip flight is a particular 
characteristic of many modern aircraft flying at high angles of attack. 
The primary causes are asymmetric forebody vortices, viscous flow 
separation of unswept and moderately swept wings, asymmetric leading- 
edge vortices of slender wings, and the strong interaction between 
aircraft aerodynamic components. 

At high angles of attack, cross flow dominates for a long pointed 
nose so that asymmetric leeward vortices or wakes will be randomly shed 
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and give undesirable side force. Even small manufacturing asymmetries 
or nose misalignment (refs. 15 and 16) could cause a cross-flow boundary 
layer to separate unevenly. Such a nonzero side force or moment occurs 
not only on the fuselage alone but also on complete airplanes, such as a 
one-tenth scale F-5E complete- airplane wind-tunnel model (ref. 17). 

Unsteady shock movement and shock boundary layer interaction could 
also cause the flow on the upper surface of the outboard section of the 
wings to separate first for aircraft with low to moderately sweptback 
wings, such as the F-4 (ref. 18), the F-5 (ref. 19), and the Gnat 
aircraft (ref. 10). Recent oil-flow studies of a single-engine general 
aviation research airplane model (ref. 20) clearly pointed out that 
asymmetric flow separation had progressed from wing-root trailing-edge 
forward and outboard of the rectangular wing. This resulted in the 
earlier stall of one wing compared to the other, consequently producing 
fluctuating pressure changes and yielding a rolling moment. 

At high angles of attack, aircraft with highly swept wings can 
profit from vortex suction lift produced by a pair of strong vortex 
cores along the upper surface of the wing, if vortex breakdown or 
vortex-sheet contact does not happen on the wing. Once the vortex 
breaks down, the vortex pattern is very sensitive to sideslip, resulting 
in a large change of the rolling moment due to a small sideslip 
disturbance. There is yet another type of asymmetric vortex pattern 
without vortex bursting. Vortex sheets from both leading edges of a 
slender wing may contact each other at zero sideslip (fig. 6) so as to 
expel reattached flow out of the center area on the upper surface of a 
wing. Possibly because of hydrodynamic (inviscid) instability (ref. 21) 
in the vortex flow field, the initially symmetric vortices are trans- 
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formed into asymmetric vortex patterns such that one vortex slides up 
over the other to produce asymmetric loads. As sweep angles of delta 
wings increased from 78 to 84 degrees, induced rolling moment at zero 
sideslip was shown (ref. 22) to rise from near zero to a very high 
value. Delta wings with leading-edge sweep angles from 80 to 86.5 
degrees (refs. 8, 23, and 24) also experienced some asymmetric vortex 
structures. Hence, delta wings with sweep angles greater than 78 
degrees ere prone to having an asymmetric flow field at zero sideslip. 

Forebody vortices and shedding vortices from upstream high-lift 
devices such as strakes and canards will heavily affect the flow over 
the wing and tail. Nonzero side forces or moments due to the strong 
interaction between aircraft aerodynamic components are evident on the 
F-5 (refs. 11 and 14) and other configurations. Analysis of wing rock 
for these configurations should account for the whole aircraft rather 
than the wing alone. 

2.2. Transient Development 

Once asymmetric flow starts, a roll-oscillation amplitude will keep 
building up if the roll damping is weak or negative. For example, an 
F-4 encountered severe wing rock in flight (ref. 25). Wind-tunnel 
force-oscillation tests (ref. 26) later confirmed that the negative roll 
damping at small amplitude near the stall angle of attack was 
responsible for transient wing rock. Other configurations such as the 
Gnat aircraft (ref. 6), the F-5 (refs. 5 and 14), the F-14 (ref. 27), a 
recently proposed fighter configuration (ref. 2), a single-engine 
general aviation research aircraft (ref. 20), the X-29A forward-swept 
wing fighter (ref. 28), a series of flying wings with vertical tails 
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(ref. 29), and a delta wing of 80-degree sweep (ref. 8), all attributed 
wing rock development to the destablizing roll damping at high angles of 
attack. For the delta wing with 80-degree sweep, the asymmetric vortex 
pattern (fig. 4) became more pronounced when the wing rolled. The 
windward vortex appeared to diffuse and shift inboard while the leeward 
vortex remained strong and moved outboard, resulting in a propelling 
rolling moment at a small roll angle. The transient wing rock due to 
roll instability at low roll amplitude will grow gradually over some 
oscillation cycles. The decisive factor in developing transient wing 
rock is the loss of roll damping at low oscillation amplitude where the 
dihedral effect is almost powerless to counteract. 

2.3. Sustaining of Limit Cycles 

At high angles of attack, the magnitude and sign of roll damping 
derivatives may vary with oscillation amplitude, frequency and roll 
rate. Several sets of wind-tunnel data, including the F-4 (ref. 26), 
the F-5 (ref. 5), a recently proposed fighter configuration (ref. 2), 
and a delta wing with 80-degree sweep (ref. 8), showed that effective 
roll damping was negative at small roll amplitude but positive at large 
roll amplitude (fig. 7). For a delta wing with 80-degree sweep, one 
reason was that an asymmetric vortex pattern (fig. 4) at large sideslip 
during rolling would yield positive roll damping, because the leeward 
vortex was displaced above and finally off the wing. Force-oscillation 
tests of this delta wing showed only little dependency of roll damping 
on oscillation frequency. However, they were conducted at frequencies 
well below wing-rock frequencies observed in free-to-roll tests 
(fig. 8). On the other hand, the roll damping of the F-4 (refs. 26 
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and 30) deteriorated quickly as the oscillation frequency was reduced. 
Cubic variation of the rolling moment with roll rate seemed to correlate 
well with wing rock for the Gnat trainer (refs. 6 and 10). This implied 
that the roll damping was a function of roll rate. Other data showing 
nonlinear dependency of the rolling moment on roll rate in rotary tests 
(refs. 2, 3, and 8) also implied that the roll damping was a function of 
roll rate (fig. 9). Thus, nonlinear roll damping depending on ampli- 
tude, frequency, and roll rate probably plays the leading role in 
sustaining the limit cycle of wing rock. 

From the viewpoint of basic fluid mechanics, the nonlinear roll 
damping probably results from various vortex asymmetries and time lags 
associated with the time history of fluid and induced mainly by the 
fluid viscosity. Convective lag is well known to exist because shedding 
vortices convect downstream at a finite speed equal to the local free- 
stream velocity and disturbances created at any point propagate at the 
local speed of sound. Because the vortex flow over a wing can not reach 
its steady-state strength and position at once in plunging motion (ref. 
31), the so-called vortex lag is created. During a plunging motion, 
pressure changes on a wing for increasing angles of attack are not 
simply reversed when angles of attack are decreasing. It was 
theoretically demonstrated (ref. 32) that vortex lag would depend on 
motion frequency and create a phase lag between the wing motion and its 
leading-edge vortex strength buildup. Development of a turbulent 
boundary layer and wake must account for the infuence of upstream flow 
history on turbulent stresses (ref. 33). Since a finite time is 
required to establish or destroy a healthy flow, a time lag in 
separation or reattachment was recorded during a boundary-layer wind- 
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tunnel test (ref. 34). This test revealed that time lag seemed to be 
larger in reattachment than in separation under conditions of cyclically 
varying pressure gradient. 

Though the roll damping is negative or weak over small oscillation 
amplitudes, it is positive at larger amplitudes for a sustained wing 
rock. Both the effective dihedral effect and positive roll damping 
(fig. 7) via aerodynamic nonlinearities at large roll amplitudes will 
gradually reduce the roll rate. As restoring moments become stronger, 
the aircraft will reach a threshold amplitude and finally reverse the 
rolling direction. Thus, wing rock is constrained to a finite amplitude 
through nonlinear roll damping. 
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3. Theoretical Development 

The main difficulty in describing the relationship between 
instantaneous aerodynamic reactions and dynamic motion variables of an 
aircraft lies in the fact that this relationship is not solely 
determined by instantaneous values of motion variables but also 
determined by the time history of the motion up to the instant in 
question. A straightforward approach is to solve the unsteady flow- 
field equations simultaneously with the equations of motion of the 
aircraft. However, to avoid the need for solving both dynamic and 
aerodynamic equations at the same time, a proper aerodynamic model 
relating aerodynamic reactions and motion variables must be developed. 

Currently, there are several aerodynamic models proposed to 
describe wing rocking. Besides some disadvantages, each of them is 
verified only for a particular configuration. 

First, it was shown that by including a cubic term in the roll 
damping derivative of a Gnat aircraft (ref. 6), the angle of sideslip 
during wing rock could be fairly well predicted by reducing the three- 
DOF lateral-directional equations of motion to a fourth-order differ- 
ential equation in 8 and by obtaining solutions from the application of 
the Beecham-Titchener method (ref. 1). When similar cubic terms were 
introduced in the static yawing moment derivative due to sideslip for 
the HP 115 research aircraft (refs. 12 and 35), the estimated amplitudes 
of limit cycles were about 40 percent greater than those obtained from 
six-DOF nonlinear simulations. However, using the nonlinear stiffness 
alone in a one-DOF model can not explain the existence of wing rock. 
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Second, nonzero static lateral-directional aerodynamic forces and 
moments at zero sideslip angle, together with aerodynamic hysteresis 
were used to reproduce the wing rock of an F-5E in a six-DOF digital 
computer simulation (ref. 11). The aerodynamic hysteresis of the 
rolling moment with sideslip was also assumed in a simple two-DOF model 
including rolling and yawing moments to demonstrate that roll hysteresis 
could be a potential cause of wing rock (ref. 36). Mathematically, such 
hysteresis (fig. 10) can be accounted for by including a function which 
can assign two possible values in the rolling moment for a given 
sideslip angle. A single-DOF model with a hysteresis loop is possible 
to yield the motion of wing rock. However, limit cycles are obtained 
only when an external disturbance is large enough to induce a sideslip 
angle to lie outside of the 0-range in the hysteresis loop. Further- 
ore, the loop shape must be determined before solving the equation of 
motion. 

Third, the variation of roll damping with sideslip angle (fig. 11a) 
such that damping was negative at small sideslip angles but positive at 
large sideslip angles could induce wing rock for a delta wing with 80- 
degree sweep (ref. 8). A one-DOF nonlinear simulation for this delta 
wing produced the motion of wing rock in close agreement with test 
results of a free-to-roll model. However, it was also shown that an 
approximate analytical solution of the one-DOF equation based on the 
Beecham-Titchener method by assuming a linear variation of roll damping 
with sideslip angle underestimated the amplitudes of wing rock by 15 
percent for angles of attack between 25 and 35 degrees. Figure lib 
indicated that in contrast to no hysteresis for static-f orce-test data, 
dynamic hysteresis was recorded during free-to-roll tests. 
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Fourth, by representing the time history effect by a lumped time 
lag (ref. 37), limit-cycle amplitudes were predicted by use of 
experimental static data (fig. 12). However, this one-DOF analysis 
could only predict the limit-cycle amplitude for delta wings with 
leading-edge sweep larger than 74°. Besides, the frequency must be 
known in advance. 

Finally, the unsteady incompressible inviscid flow equations and 
the one-DOF equation of motion were simultaneously solved in reference 
38. It resulted in dynamic hysteresis of rolling moment versus bank 
angles (fig. 13). However, it failed to predict the maximum limit-cycle 
amplitude for an 80-degree delta wing because the vortex breakdown 
effects were not taken into account. Further, this method could not 
provide any aerodynamic derivatives causing wing rock. In addition to 
some numerical simulation problems, it is very time-consuming. 

To describe the wing-rock characteristics of any aircraft, the 
aerodynamic model developed here will be more general to include all the 
essential aerodynamic nonlinearities presented in the aforementioned 
current aerodynamic models. In the meantime, the nonlinear roll damping 
depending on amplitude, frequency, and roll rate may be a primary factor 
in sustaining wing rock and should be accounted for properly. 

Since wing rock is dominated by the roll oscillation, three-DOF 
lateral-directional equations of motion will be used in the present 
investigation, instead of the general six-DOF equations of motion. 

Accurate solutions can be achieved through numerical integrations. 
However, purely numerical solutions do not provide direct physical 
insight into the dynamic mechanisms of wing rock. A useful analytic 
approximation to the characteristic motion of wing rock is more suitable 
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for a parametric study of aerodynamic and dynamic causes. The Beecham- 
Titchener method (ref. 1) will be applied and extended here for this 
purpose. 

3.1. One-Degree-of-Freedom Model 

A one-DOF rolling equation is to be studied first to describe the 
idea. By ignoring any pitching or yawing effects, the principal motion 
variable in pure rolling motion is the Euler roll angle, ij>. Flight path 
properties such as flight velocity, altitude, Mach number, Reynolds 
number, and mass distribution are all assumed to remain constant 
throughout the motion. Effects of elastic, gravitational, and 
propulsive forces are all excluded. For a rigid aircraft with the 
observer fixed in the geometric body-axis system (fig. 14) whose origin 
is at the aircraft center of mass, the equation of motion is given by 

^xx * , , , , . 

p = C. (t) (1) 

qSb 

where I xx is the rolling moment of inertia, q is the dynamic pressure, S 

is the reference wing area, b is the wing span, t is the time variable, 

• 

P is the angular acceleration about the longitudinal axis, and C^(t) is 
the total aerodynamic rolling moment coefficient. Equation (1) simply 
states that the inertial rolling moment is balanced by the time- 
dependent aerodynamic rolling moment. Since classical stability 
derivatives in Taylor series expansions alone are not adequate to 
account for the aerodynamic nonlinearities at high angles of attack, the 
concept of a simplified nonlinear functional (ref. 39) is used to 
describe the total aerodynamic rolling moment. 
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After examining the wind-tunnel data of a flat-plate delta wing 
with 80-degree sweep (ref. 8) subjected to static-force, forced- 
oscillation, rotary, and free-to-roll tests, the characteristic 
variables of C (t) are selected to be a (the steady-state angle of 

& S 

attack), g (sideslip angle), p (the reduced rolling velocity about the 
X-axis), and 6^ (aileron deflection). Hence, the total aerodynamic 
rolling moment coefficient at fixed a g can be expressed in the form of a 
functional as 


C (t) = C £ (a g , 8(t), p(t), 6 £t)) (2) 

Practically, the functional will be replaced by a multivariable 
function which depends on a few characteristic parameters rather than on 
continuous functions. Furthermore, the multivariable function will be 
expanded into a form of component buildup as 


where 


c e (t) ' + V + V 55 * V 6 a 

0 8 pt S A 


c, - c + C 1 3 1 + c [ p I 

pt X p0 x p8 PP 


(3) 

(4) 


The total aerodynamic rolling moment coefficient in equation (3) is 

composed of (1) C 5 (a ,8=0), the aerodynamic rolling moment 
*0 3 

coefficient arising from asymmetric flow at zero sideslip angle, being 

highly dependent on aircraft configurations; (2) Q (a , 8), the 

p 

dihedral effect, being usually negative (statically stable), but 

possibly positive (statically unstable) for some configurations (refs. 

40 and 41) at high angles of attack; (3) C. (a , 8, p) , the total or 

Xj s 
pt 
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effective aerodynamic roll damping coefficient, being a strong nonlinear 
function of a g , 3, and p; and (4) (a g , P> <$ ) > the aerodynamic 

5 a 

roll derivative due to aileron deflection. 

Equation (4) represents the present proposed aerodynamic model to 

describe the aforementioned aerodynamic nonlinearities. This expression 

is expected to give negative roll damping for small roll amplitudes and 

positive roll damping for large amplitudes in accord with experimental 

data. If is plotted against 3 or p, the resulting curve is assumed 
P 

to be symmetric with respect to 3 = 0 or p = 0 respectively. The 
absolute signs in equation (4) will cause the values of C at 
different instants to be the same if |3| or |pj have the same values at 
the corresponding instants in a limit-cycle oscillation. At a 
fixed 3, may have two values if the sign of p is different. 
Furthermore, at a fixed p, may also have two values if the sign 
of 3 Is different. Thus, 3- and p-dependent roll damping is able to 
produce discontinuous and double-value rolling moments with respect 
to 3 or p. In other words, it may result in a hysteresis loop of 
rolling moment with respect to 3 due to the cyclical change of total 
roll damping. 

Equation (1) together with the kinematic relationships 

p = <J> (5) 

3 = sin ^ (sin a sin 4) (6) 

s 

form a nonlinear second-order differential equiation in <j>, with 
coefficients being also functions of time. To simplify equation (3), 
equation (6) is approximated by 
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( 7 ) 


3 = d> sin a 
r s 


This is a good approximation, because 3 of equation (7) is about 1 
degree higher than 3 of equation (6) for most cases of wing rock. 
Substituting equations (3), (4), (5), and (7) into (1) yields 


I 

xx 

qSb 




C. <5. + sin a C. | + 
£ A s 

5 A 6 


+ 



+ sin a 

s 



M + c 4 

pp 


1 4>b | j>b 
' 2V 1 ' 2V 


( 8 ) 


For convenience, equation (8) is rewritten in terms of dimensional 
derivatives as 


= L_+ L. 6 + sin a L fl $ + 

0 6, A s 3 

A 


• • 


+ (L p0 + sin a g L pg |*| + L pp |<fr|) * 


(9) 


where dimensional derivatives are defined by: 


q S b _ , -2., 

L - V - C. («<= > 

xx 0 


s - v (sec ' 2) 

xx 3 


T q S b^ „ , -1, 

L P 3 = 2VI C £ . (S6C 5 
v xx p3 


PP 


I 

XX 

* 

- q s 

b 2 

2 V 

*xx 

to 
1 O- 

II 

b 3 

4 V 2 

X xx 


- 2 , 


C. (sec *) 


(10) 


Z % (dimensionless) 
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Equation (9) is the nonlinear equation of motion to be solved for the 
one-DOF model. 

3.2. Approximate Analytic Solutions of the One-Degree-of-Freedom Model 

To solve equation (9) so that the amplitudes and periods of wing 
rock can be determined, the Beecham-Titchener method (ref. 1) appears to 
be the most appropriate. This method is basically an extension of the 
classical Krylov-Bogoliubov method. It is shown to be a combination of 
the averaging principle and the method of variation of parameters (ref. 
42). For a damped linear system (ref. 1), the first-order approximation 
of the Beecham-Titchener method gives the exact solution. For such a 
system, the Krylov-Bogoliubov method yields an incorrect undamped 
natural frequency. Other examples (ref. 1) involving nonlinear terms 
are solved with excellent accuracy by using the Beecham-Titchener 
method. This method is also capable of obtaining the higher-order 
approximations. Therefore, this method is applied here to obtain an 
approximate analytic solution to equation (9). 

By analogy with the solution for a linear second-order differential 
equation, the time-dependent solution of equation (9) is assumed to be 

<j>(t) = a(t) cos v(t) (11) 

That is, a product of an amplitude function a(t) and a sinusoidal 
function in terms of an argument v(t) is assumed. Both u> and X are 
defined in such a way that 
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Ui 


V 


( 12 ) 


X = 


a/a 


(13) 


To the first-order approximation (appendix A), equation (9) without 
control deflection is reduced to 


X^-w^ = sin a L„ + XL „ + aX sin a L „ + -f— a X u> L 

s 8 pO 3tt s pB 3tt pp 


(14) 


2X(u = to L _ + ~ a (o sin a L . + ~ a (X^ + 2w^) L 

pO 3ir s pp 3 tt pp 


(15) 


Equation (9), a formidable nonlinear second-order differential equation, 
is thus reduced to two coupled nonlinear first-order differential 
equations. Both X and w turn out to be functions of the amplitude. 
Equations (14) and (15) can be further integrated by a numerical 
technique. All the aerodynamic derivatives in equations (14) and (15) 
should be evaluated at some average conditions according to the Beecham- 
Titchener method. 

If 52, P, and A are the circular frequency, the period, and the 
amplitude respectively for a steady-state limit-cycle wing rock, then 
they can be obtained by setting X to zero in equations (14) and (15). 

The closed-form solutions can be shown (appendix A) to be 


52 = (-sin a 

s 


V 


1/2 


/ 9 S b . „ ,1/2 

= (- :1 - sin a C ) 

I s Z„ 

xx B 


(16) 
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p 


(17) 


2tt 

a 


a = - 


3n jpO 

sin a L . + 2 0 L 
s pg pp 


3tt 

4 


'pO 


sin a s C * + ~ C i 

pg PP 


They must be real and positive to be physically realistic. As can 
seen from equations (16) and (17), other than inertial properties, 

limit-cycle period depends only on C . The limit-cycle amplitude 

8 

function of C. , C , and c. . Equation (18) can be reduced to 
pO pg pp 

expression used in reference 8 if C is ignored. 

PP 


(18) 


be 
the 
is a 
the 


3.3. Three-Degree-of-Freedom Model 

To obtain a better understanding of wing rock, a three-DOF model 
will be developed. Since the angle of attack for a sustained wing rock 
remains essentially constant, the effects of the longitudinal mode on 
lateral-directional modes can be neglected in the theoretical analysis. 

Most aircraft are symmetric with reference to a vertical plane (X-Z 
plane) aligned with the longitudinal axis (X-axis) so that I xy = 0 and 
Iy z = 0 are satisfied. If inertial coupling moments associated with pq 
or qr are ignored, the equations of motion for three-DOF lateral- 
directional modes (ref. 43), based on body-fixed axes, are 


(v - wp + ur - g sin <(> cos 0) = Cy(t) 
qS 


*xx * ^XZ ‘ „ , S 

— p _ __ r = c £ (t) 

qSb qSb 


(19) 


( 20 ) 


19 



( 21 ) 




" C n (t) 


where m is the airplane mass and g is the acceleration of gravity; u, v, 
and w are velocity components along X, Y, and Z axis; p and r are roll 
and yaw rates about X and Z axes respectively; 9 is the pitch attitude 
angle; Cy and are the total aerodynamic side-force and yawing moment 
coefficients respectively. C Y ( t )» C £ (t), and C (t) can be written in 
forms of functional as 


C y (t) - Cy (a g , 8(t) , p(t), r(t), 

C £ (t) = C £ (a g , 6(t) , p(t), r( t) , 

C (t) = C (a , 6(t), p(t), r(t) , 

n ns 


« A <t>, 

V‘» 

(22) 

« A (0, 

v*» 

(23) 

6 A (t), 

V'” 

(24) 


where r is the reduced yaw rate and <5^ is the rudder deflection angle. 
Similar to the one-DOF analysis, Cy(t), C^(t), and (^(t) are further 
expanded into forms of component buildup as 


C Y (t) " V + C Y 8 + S P + C Y ~ V + C Y . 6 A + V 6 R (25) 


0 8 


P r 5 6 

y A R 


C £ (t) = C f n + V 3 +( S n + J B ! +C £ I ? I )P + 

06 pO p8 pp 


C * r + V S A + V S R 
r 5 A S R 


(26) 


c „ (t> * V + V 6 + < C „ „ + c n l«l + C n I5l> P + 

0 6 pO p$ pp 


C r + C 6 a + C 6_ 
n n. A n. R 

r o » 

A R 


(27) 
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Because the motion of wing rock is primarily rolling around the longi- 
tudinal body axis, the total aerodynamic yawing moment due to roll rate 
in equation (27) is expected to be as important as the total aerodynamic 
rolling moment due to roll rate in equation (26). Both roll and yaw 
damping terms associated with sideslip and roll rate are thus assumed to 
have similar forms. 

To simplify the equations of motion, the following approximations 
are applied. 


sin <j> » <J> 

COS (j) s> 1 

u/V a cos a (28) 

s 

v/V « B 
w/V « sin a 

s 


Equations (25) - (27) are now rewritten as 


a - 


p sin a + r cos a - §■ <P cos 9 = ■ a — (C„ + C„ 3 + C p + C r 

s s V mV Y_ Y„ Y Y 

0 B p r 


_ qs 


+ V «A +< V. V 

5 a 5 r 


(29) 


p _ JE5. r = SSb [ C + c £ B + (C z + C z |e|+C £ |p|)p + C ? 

xx 0 B pO pB pp r 


xx 


+ C 6 + C. 6 1 

£. A i. R J 

6 a 5 r 


(30) 


r - 


I s ? - f 1 < C » + C n 0 6 + <C n „ + C n J 6 ' + C „ 151’? + C „ ? 

zz zz 0 8 pO p8 PP r 


+ V 4 A + V V 

4 A 4 R 


( 31 ) 
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Rewriting equations (29), (30), and (31) in terms of dimensional 
derivatives yields 


B - 


p sin ct g + r cos « g - | <f> cos 0 = Yq+Y^+Y p+Y r r*Y fi fi A +Y fi 6 R (32) 


A “ R 


P - 


— r ■ l 0 + V + a p0 +L p8 |S|+l pp!' l|) p + irl + h S A + h 5 R <33) 

xx cosa A R 

s 


r - 


T 2 - p - VV +< VV B| %> i)P+ -Sr™6 V H « S R < 34 > 


Z2 


cosa A R 
s 


where 


Y = c (sec”" X ) 
0 mV L Y q ^ seC ; 


Y =s -SA q (sec~^) 
8 mV Y0 U C ; 


Y « — - b ... C (dimensionless) 
P 2mV^ Y p 


Y = -S ^ iL. c (dimensionless) 
r 2raV 2 Y r 


Y S. - 2v V <3e '"‘ ) 

4 S A 


T «„ ’ « V (se =' I > 

R S R 


L o ” i 55 - V (sec * 2) 

xx 0 


L s ■ c t (sec ' 2) 

xx g 


T qSb 2 _ . —I, 

L p0 21 V C fc (seC } 

xx pO 


L p 6 ’ 2T-V C t , (sec '‘ ) 05) 

XX p8 


qSb^ 

L = —3- T C. (dimensionless) 

pp . T ,,2 Z 

41 V pp 
xx 


qSb 2 , -1 N 

L r = 2T~V C £ C0Sa s’ (sec > 
xx r 


<sep ' 2) 

A 


L fi = C (sec 2 ) 
S R X xx \ 


N 0 = C n ^ SeC ^ 
0 I n 

zz 0 


N = C (sec 2 ) 

8 : p 2 n s 
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qSb 2 

21 


zz 


V 



(sec *) 



qSb 2 

21 


zz 


V 



(sec 1) 


N = —3 T C (dimensionless) N = -i= — - C cos a (sec ) 

pp , _ „2 n r 21 V n s 

v 41 V pp zz r 

zz 


qSb 


zz 


(sec ) 



(sec ) 


The equations of motion for the three-DOF model are finally obtained by 
rearranging equations (32) - (34) as 


M Y. + YJ + Y p + Y r + Y * + Y. 5, + Y. 

0 3 p r <f> T 6 A 6 R 

A 


(36) 


p - E 0 + V + <%0 + Sb | b I + Spl»l) » + iii- + V 4 + V R 

S 


(37) 


* - 5 0 * V + <S pO + v' 61 + S pptpl> p +£&+ V S A + S « R (38) 


where 


Y = Y 
0 0 


Y, = Y 

3 6 


Y = Y + sin a 
p p s 


Y =(Y - cos a )cosa 

r v r s s 


\ - t c ° 3 8 


Y = Y 
6 5 

A A 


Y = Y 

6 R 6 r 


i = a - y~ X t > 1 

xxzz 
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3.4. Approximate Analytic Solutions of the Three-Degree-of-Freedoro Model 
The Beecham-Titchener method was originally developed for a single 
nonlinear second-order differential equation only. This method will be 
extended here for the three-DOF equations of motion to obtain the 
approximate analytic solutions for wing rock. The solutions to 
equations (36) - (38) are assumed to be 

0(t) = a(t) £ (t) cos [ v ( t) + e(t)] (42) 

<j)(t) = a(t) cos v(t) (43) 

if<(t) = a(t) n(t) cos [ v ( t ) + 5 ( t ) ] (44) 


where £ and n are the instantaneous amplitude ratios of 0 and ^ modes to 
the <(> mode and e and 6 are the phase angles of 0 and modes with 

respect to the <(> mode. 

To the first-order approximation (appendix B) , equations (36) - 
(38) without control deflections can be shown to be 


C(Xcose - uisine) = Y„ £cose + Y X + Y n (Xcos6 - wsind) + Y, (45) 

0 p r (p 


C(Xsine + rncose) = Y £sine + Y u> + Y n (Xsind + wcosd) (46) 

p p r 


2 2- — g _ g — 

X -m = L Ccose + L „X + — L „ a£X + — L aXrn 

0 pO 3 tt p0 3 it pp 


+ L r n (Xcosd - oisind) 


(47) 
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( 48 ) 


2Au) = L a ^sine + L -co + 4 — L . a ?i») + 4 — L a(A^ + 2a>^) 

8 pO 3 tt p8 3 it pp 

+ L r n (Asin6 + ujcosS) 

n t (A 2 -<j 0 2 ) cos 5 - 2Au)sin6] = N. Ccose + N A + 4— N „ a£A 

8 pO in p8 

+ 4— N aAto + N n(Acos6 - msin6) (49) 

Jr pp r 

2 2 - - 4 - 

n [ ( A - 0 ) ) sin6 + 2 Aujcos 6] = N. 5sine + N „u> + -x— N . a£iu 

8 p(J JTT p8 

+ 4— N a(A^ + 2u>^) + N n (Asin6 + ucosfi) (50) 

Jit pp r 

Again, all the aerodynamic derivatives in equations (45) - (50) should 
be evaluated at some average conditions according to the Beecham- 
Tltchener method. These six coupled nonlinear first-order differential 
equations can be integrated numerically to find a, £ , n , e, v,and 5 as 
functions of t for known initial conditions. 

To suppress the wing rock effectively, the main aerodynamic 
nonlinearities causing the self-induced limit cycle must be understood 
and justified first. An approximate theory, such as the one developed 
here, is an efficient tool for recognizing the aerodynamic nonlineari- 
ties important to the amplitude and frequency of a limit cycle. The 
amplitude of a limit cycle is the primary concern as a handling-quality 
or safety problem. Therefore, the main effort here is to find out the 
effects of amplitudes of 8 and ^ modes on that of the <j> mode for a 
steady-state limit-cycle oscillation. 
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To obtain the limit-cycle amplitude of a sustained wing rock, X is 
set to zero. In addition, the phase angles, e and 6, are ignored since 
they are found numerically to be insignificant for a steady-state 
constant-amplitude wing rock. Equations (45) - (50) are then reduced to 


0 = V + " Y * 


G = Y + Y H 
P r 


a - V 


0 = L + ^r— L Q AG + -i — L Afi + LH 
pO 3ir p8 3ir pp r 


HIT = -N G 

p 


- 4 - 8 - 

0 = N n — N . AG + -t— N A!i + N H 
pO 3ir pB 3 tt pp r 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


where G and H are the amplitude ratios of 8 and <|> modes to the 4> mode at 
limit-cycle conditions. 

It is assumed that the frequencies of <f> , B,and i|/ modes are the same 
in equations (45) - (50). Because equations (51) and (55) are related 
to the frequencies of 8 and iJj modes respectively, they can be ignored in 
order to find the effects of 8-ty amplitudes on the <f> amplitude for a 
steady-state wing rock. The limit-cycle conditions are then obtained in 
appendix B by rearranging equations (52), (53), (54), and (56) into the 
forms 
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G = Y + HY 
P r 


= sin a - H cos^a + qS ^ -~ C + H C cos a 

*2 Y_ '>_*tt 2 Y_ s 


(59) 
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The roll amplitude of a limit cycle is affected not only by C , 

C. , C„ , and C as in the one-DOF model but also by c , Cy , Cy , 
pO p3 pp n 6 P r 

C , C ,C , C„ , and C as shown in equation (57). Equations 
pO pg pp r r 

(57) - (60) are coupled nonlinear algebraic equations in terms of A, G, 
H, and ft . They can be solved through simple iterations started with G 


= sin a g . For a steady-state wing rock, all solutions must be real and 
positive. 


* 
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4 


Numerical Results and Discussions 


It would be of great interest if the predictions of the present 
theoretical results could be shown to agree with measured wing-rock 
characteristics of aircraft, such as the F-4, F-5, F-14, X-29A, Gnat, 
Harrier, HP 115, and some reentry vehicles. Unfortunately, not all 
aerodynamic derivatives needed in equations (57) - (60) are available. 

In addition, it is difficult to estimate the high-a aerodynamic 
derivatives for a full-scale airplane especially when effects of vortex 
breakdown, viscous flow separation, and the aerodynamic hysteresis on 
the wing become significant. 

A recent theoretical study of slender wings (ref. 44) based on the 
quasi-vortex-lattice method (ref. 45) and the suction analogy (ref. 46) 
showed that the longitudinal and lateral-directional aerodynamics with 
vortex-breakdown effects at high a can be reasonably well predicted. 
Therefore, the updated steady-flow aerodynamics computer code 
(VORSTAB)** based on reference 44 will be used here to calculate all 
required aerodynamic derivatives for slender wings for the present 
purpose. 

To test the one-DOF theory, a delta wing with 80° leading-edge 
sweep (ref. 8) is chosen. To verify the three-DOF theory, eight flying 
wings with vertical tails (ref. 29) are selected. 

Having formulated the mathematical model on body-fixed axes, 

6 derivatives must be included in all p and r derivatives. However, no 
existing theoretical methods are available for predicting 8 derivatives. 

**The VORSTAB computer code is available through the NASA Langley 
Research Center. 
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Since the aerodynamic nonlinearities of wing rock mainly arise from the 

nonlinear roll damping of the wing, it is decided to correlate the roll 

damping derivatives (C^ ) based on stability axes and calculated from 

P s 

the VORSTAB code with the roll damping derivatives (C + C . sin a ) 

jC /C _ s 

P 3 

based on body axes and obtained in forced-oscillation tests. The reason 

is that the effect of C. . on C 5 based on body axes is similar to that 

3 P 

of C. on C based on stability axes. The results of four delta 
X* X* 

r p 

s s 

wings with 80°, 76°, 74°, and 60° leading-edge sweep (refs. 8, 22, 

47,48, and 49) are presented in figures 15-18, Two conclusions can be 

made with regard to comparison between the calculated 

C. and data of C C . sin a . First, the agreement is reasonably 
Xf x» • x> n s 

p s P 3 

good if the amplitude in forced-oscillation test is not too small and 
the vortex-breakdown effect is not significant. Second, the agreement 
is better if data of the forced-oscillation technique are compared 
instead of those of the free-to-damp-oscillation method. 

To simulate a one-DOF free-to-roll test or forced-oscillation test 
based on body axes, an appropriate location of moment center must be 
chosen if the VORSTAB code is to be used to calculate the aerodynamic 
derivatives on stability axes. It should be noted that the VORSTAB code 
is to simulate the rotary-test results. The moment-center locations 
used in the present study are set at the aerodynamic-center location 
taken from wind-tunnel data (refs. 22 and 47). Figures 15-18 confirm 
that the aerodynamic-center location is a good choice for the moment 


31 



center if the aerodynamic derivatives are to be calculated on stability 


axes. Furthermore, numerical results also show that the aerodynamic 
derivatives change very little if the moment-center location moves +0.10 
c from the aerodynamic-center location. 

Knowing that the roll-dominated wing rock is mainly a lateral 
problem caused by the wing, the VORSTAB code is now used to calculate 


and for (actually c +C^. sin a s ) and 

p s p p 8 

C 0 (a tually C. -C . cos a ) respectively. The other aerodynamic 
r r n 0 S 

derivatives are calculated on body axes. 


4.1. One-Degree-of-Freedom solutions 

To verify the one-DOF theory, a flat-plate delta wing with 80° 
leading-edge sweep is chosen because its wind-tunnel data of static- 
force, forced-oscillation, rotary, and free-to-roll tests have been well 
conducted and published (ref. 8). 

According to the Beecham-Titchener method, all aerodynamic 
derivatives should be evaluated at some average dynamic conditions. 

Since the interaction between aerodynamics and dynamics in wing rock is 
strongly nonlinear, an iterative scheme (appendix C) is developed to 
calculate the average aerodynamic derivatives and dynamic characteris- 
tics for a steady-state self-induced limit cycle. In other words, the 
aerodynamic derivatives are calculated at some average dynamic 
conditions (i.e. certain average values of 8 and p) which in turn depend 
on the aerodynamic characteristics. 
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The calculated longitudinal aerodynamic derivatives of the test 

configuration are presented in figure 19. The predicted results agree 

quite well with wind-tunnel data even at a beyond a (a = 38°). 

The limit-cycle amplitude is the main concern in studying the wing- 

rock problem. The calculated results are plotted in figure 20. To the 

first-order approximation, the prediction is generally good up to a of 

38°. Beyond that a, a converged solution is difficult to obtain. 

Besides the strong vortex breakdown effects at higher a, the viscous 

flow separation will probably also become important. The latter effect 

is not included in the VORSTAB code. The predicted limit-cycle period 

is reasonably good as shown in figure 21. 

For a limit-cycle oscillation, the values of C are positive and 

P0 

those of C are negative; while those of r may be either positive 

P0 PP 

or negative (fig. 22). The effect of C is to decrease the predicted 

PP 

limit-cycle amplitudes for a below 30° and above 36°; while it increases 

the predicted amplitudes for a in between. Without C , the predicted 

PP 

amplitude is 54.4° (overestimated) at a of 25° and is 31.1° (under- 
estimated) at a of 32° . The predicted maximum amplitude with C is 

Xr 

P p 

40.7° at a of 32° and is very reasonable when compared with available 
data (fig. 20). Thus, the roll-rate dependency of roll damping plays an 
important role in determining the limit-cycle amplitude. 

The theoretical variations of all terms in equation (4) for a = 27° 
and 32° are plotted in figures 23 and 24 respectively. It is seen that 
magnitudes of and |p| ar® much higher than those of J g | at 

p0 pp p£ 

lower bank angles while C | B | is dominating at higher bank angles. 

P0 

The shape of C and its values are affected by the sign of c. 

x> . x> 

pt pp 
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The histograms of versus $ for one limit cycle of wing rock 
based on equations (3) and (4) without C or the control term at a of 
27° and 32° are presented in figures 25 and 26 respectively. The loops 
give evidence of dynamic hysteresis during the wing rocking. Agreement 
of theoretical contours with the data derived from free-to-roll tests 
(ref. 8) is reasonable. The theoretical loops are symmetric against 
positive and negative bank angles, but the experimental loops depict 
some asymmetries. 

To understand the loop asymmetries, it is known that the induced 

rolling moment at zero sideslip C increases with increasing a for a 

l 0 

highly swept delta wing. The existing data of C (refs. 8, 9, and 22) 

£ 0 

vary with the model and the wind tunnel used in the static-force tests 

for a delta wing with 80° sweep. The vortex asymmetry is the main 

factor to cause the induced rolling moment at zero sideslip. The net 

difference of experimental C between clockwise and counterclockwise 

loops (i.e. AC^ at <j> of 0°) determines which direction the wing will 

tilt and how the dynamic loop will shift with respect to the theoretical 

prediction. The VORSTAB code can not predict C. . At a of 27°, AC, is 

0 0 

+0.01 (fig. 25). The dynamic loops shift up. The positive amplitude is 

5° greater than the negative amplitude. At a of 32°, AC is -0.05 

*0 

(fig. 26). The dynamic loops rotate counter-clockwise. The positive 
amplitude is 5° less than the negative amplitude. In addition, during 
the oscillation the effect of time lags strongly affects the vortex 
patterns, causing the dynamic hysteresis. It is the combination of 
vortex asymmetry and time lag which produces the asymmetrical histogram 
in experiment. 
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The total aerodynamic energy exchanged over a limit-cycle 
oscillation is given by 


AE = isb / C (t) 4> ( t ) dt (61) 

C(t) 
or 

AE = qSb / C (<}>) d(j> (62) 

C(<j>) * 


where C(t) and C(4> ) are the contours obtained by plotting C versus t 
and <|> respectively. For a steady-state wing rock, C(<j>) is a closed 
curve over one oscillation cycle as shown in figures 25 and 26, where 
the arrows indicate the direction of increasing f'me. As derived in 
appendix D, equation (62) can be reduced to 


AE = 2qSb 2 ft/V / C (A 2 - <J> 2 ) 1/2 d<j> (63) 

0 pt 


This equation points out that the net energy transferred is 

exclusively determined by the total aerodynamic roll damping. In other 

words, to the first-order approximation it has nothing to do with 

C and C . The areas within the clockwise loops are positive or 
£ 0 l S> 

dynamically undamped due to a positive C so that energy is added to 

Pt 

the aircraft. On the other hand, the areas within the counterclockwise 

loops are negative or dynamically damped due to a negative C so that 

Pt 

energy is consumed by the aircraft. It is shown in appendix D that AE 
is theoretically zero, or the area of the destabilizing loops equals 
that of the stabilizing loops. Therefore, the net aerodynamic energy 
added to the aircraft and extracted from it is theoretically zero over 
an ideal limit-cycle oscillation. Numerical results (figs. 25 and 26) 
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further testify that the energy balance between the work input and the 

energy consumed is required to sustain a steady-state wing rock. 

The intersection of stabilizing and destabilizing loops represents 

a turning point of roll damping because C is zero at that point. The 

Pt 

critical bank and sideslip angles, 4> c and 0^, are defined as the angles 
of zero total aerodynamic roll damping. The aircraft is dynamically 
undamped below 4> ^ and damped above The solutions of $ are (appendix D) 


!<f> c l 3 ^7 t" s + (s 2 - 4rd) 1/2 ] 
U I = - — = - 22 

|Y c' 2r sin a C„ 

|<t> c ! = [-S - (s 2 - 4rd) 1,/2 ] 


if 


if 


if 


C i > ° 
PP 

C * = ° 
PP 

C * < 0 


PP 


(64) 


where 


, . „ .2 , ,bfl _ .2 

r - (sin n s ^ > H- (jf C, ) 

PB PP 

s = 2 sin a C C (65) 

s Z n Z . 
pO p0 

. = r 2 A r ^ 

d C Z n ~ ( 2V A C l ) 

pO pp 

The predicted 4» and 0^ are shown in figure 27. It is fairly good 
as compared with the available data (ref. 8). 
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4.2. Three-Degree-of-Freedom Solutions 

Simple flying wings with a vertical tail but without a horizontal 
tail or fuselage were extensively tested in the NACA Langley free-flight 
tunnel (ref. 29). Four models having delta wings with 53°, 63.5°, 76°, 
and 82.9° leading-edge sweep and four cropped delta wings with a taper 
ratio of 0.5 are used in verification of the three-DOF theory. 

The iterative scheme used in the one-DOF theory (appendix C) is 
also applicable here. The theoretical results of two flying wings are 
presented in figures 28 and 29. F or both configurations, a 76°-delta 
wing and a 76°-cropped-delta wing with a vertical tail, the predicted 
lift coefficients agree well with data up to a of 30°. 

The three-DOF theory overpredicts the starting a of wing rock by 
about 5° as shown in figures 28 and 29. The shaded areas in both 
figures represent the approximate a range in which the steady-state wing 
rock was observed. Both a ranges extend over 5° of a approximately. 

At a below that range, no wing rock was observed. At a above that 
range, a constant-amplitude roll oscillation kept increasing rapidly 
as a was increased further. As a result, the flying model rolled 
completely over out of control. The predicted a ranges of wing rock for 
both configurations are less than 10° of a. The solid lines in both 
figures can not be extended to higher-a values because the iterative 
scheme (appendix C) yields no converged solutions or equations (57) - 
(60) gives negative amplitudes or imaginary frequencies. Since no 
limit-cycle amplitudes had been recorded during the flying tests, it is 
impossible to verify the theoretical amplitudes for both flying models. 

Before the flying model reached the a of maximum limit-cycle 
amplitude, it had probably rolled off because of the great reduction in 
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the elevon effectiveness at high a. Furthermore, the above theoretical 
results are calculated for free oscillations without control deflection 
while the elevens had been used during the free-f light-tunnel tests. At 
any rate, additional data on these configurations are needed to 
ascertain the accuracy of the present theory. 

Only five flying models with an aspect ratio equal to or less than 
1 were reported to exhibit wing rock. The theoretical prediction (fig. 
30) has confirmed the same trend. Agreement between the predicted 
onset a of limit cycle and data is within ±5°. Since the free-flight- 
tunnel models simulate the more realistic cases as compared to one-DOF 
free-to-roll models, figure 30 can serve as a preliminary design guide 
for aircraft design engineers. 

The predicted period of a limit cycle decreases as a increases. 

The calculated limit-cycle periods are generally greater than four 

seconds for those five flying models exhibiting wing rock. No data of 

frequency or period were reported. Hence, no comparison is made. 

The amplitude ratio of $ to $ modes at the limit-cycle condition G 

varies from 0.01 to 0.30 for the above five flying models at 

different a. On the other hand, the amplitude ratio of to <j> modes at 

the limit-cycle condition H remains approximately 0.11 for most cases. 

Therefore, 8 and ip modes have some effects on the roll amplitude of 

limit cycle in the three-DOF theory. 

It is found in the above numerical analyses that in determining the 

wing-rock characteristics the importance of C, , C_, , and C is next 

r p0 r 

to that of C , C , and C ; while the effects of r , r , C , 
p0 pB PP P r p0 

and C are almost negligible. In the one-DOF theory, it is required 
"PP 

that C be positive (undamped) to cause a limit-cycle motion. In 
PO 
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contrast to the one-DOF theory, wing rock, is possible in the three-DOF 


theory even if C is negative (damped) but weak. The reason is that 
PO 

the combined effects of C. , C , and C may make the aircraft more 

r pO r 

unstable in roll oscillations. 

In general, there are three possible approaches to suppress or 
prevent wing rock: (1) airframe reshaping, (2) maneuver limiting, and 

(3) employing a stability augmentation system (SAS) or an automatic 
flight control system (AFCS). The first approach requires the detailed 
aerodynamic reshaping of the basic airframe configuration to provide 
inherently wing-rock-free capability. The second approach gives a fast 
solution eliminating wing rock by simply adopting an a-limiter. 

However, it may degrade usable maneuverability. Thanks to tremendous 
advances in avionic technology, the third approach employing SAS or AFCS 
has become the most effective method for attaining strong resistance to 
wing rock without degrading high maneuverability. 

Because wing rock is mainly due to the loss of roll damping at 
high a, it can be effectively suppressed by using a roll damper SAS to 
artificially augment roll damping. For example, a 76°-delta wing with a 
vertical tail is predicted to exhibit a limit-cycle amplitude of 50° 
at a of 20°, but the limit cycle of this flying wing can be eliminated 
completely by incorporating an ideal roll damper as shown in appendix 
E. 
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5. 


Conclusions and Recommendations 


5.1 . Conclusions 

A unified theory based on extensive experimental data is proposed 
to elucidate the basic mechanisms of wing rock. It is concluded that 
wing rock Is triggered by flow asymmetries, developed by negative or 
weakly positive roll damping, and sustained by nonlinear aerodynamic 
roll damping. 

A new and more general nonlinear aerodynamic model that includes 
all essential aerodynamic nonlinearities causing wing rock has been 
developed. The key element is the total aerodynamic roll damping which 
depends on oscillation amplitude and roll rate and is formulated to 
encompass the aerodynamic hysteresis implicitly. 

The Beecham-Titchener method has been applied and extended to solve 
the one- and three-DOF equations of motion. To the first-order 
approximation, the equations of motion are reduced to coupled nonlinear 
algebraic equations at limit-cycle conditions. These equations are then 
solved through iterations for the limit-cycle amplitude and frequency. 

On the other hand, closed-form solutions are obtained for the one-DOF 
theory. 

To determine the numerical values of the limit-cycle amplitude and 

frequency, the important aerodynamic derivatives such as C , C , C , 

%0 pB %p 

C„ , C , C , C. , and C , are required. Because of the strong 
r pO r 8 8 

interaction between aerodynamics and dynamics during wing rocking, these 
derivatives are evaluated with a steady-flow aerodynamics computer code 
at some average dynamic conditions in an iterative manner. 

Theoretical results show that to sustain a steady-state wing rock, 
the total aerodynamic roll damping must be negative at small bank angles 


40 


in the one-DOF case but it may be negative or weakly positive in the 
three-DOF case. On the other hand, at larger bank angles it must be 
positive in both cases. Furthermore, the one-DOF theory indicates that 
the energy balance between the work input and the energy consumed during 
motion is another reason for a self-induced wing rock. 

Good agreement between theoretical and experimental results has 
been obtained for many slender wings, of which all but one are equipped 
with a vertical tail. 

5.2. Recommendations 

Good comparison of theoretical results with data has been shown for 
some simple aerodynamic configurations. However, additional comparisons 
should be done for more configurations if data are available. 

More experiments should be conducted to understand some fundamental 
dynamic phenomena. First, systematic forced-oscillation and rotary 
tests should be made to investigate the effects of amplitude, frequency, 
and roll rate on the roll damping for simple configurations at high 
angles of attack. Second, curved-flow and rolling-flow and pure 
lateral-translation-oscillation techniques should be applied to deter- 
mining the pure derivatives such as C , C , C , C , C . , and C . . 

x x n n X— n . 

p r p r 0 0 

Finally, since the one-DOF free-to-roll test seems inadequate to predict 
wing-rock characteristics in free flight, more free-f light tests should 
be conducted. 

At present, it is very difficult to predict theoretically the 
high-a aerodynamic derivatives for a complete airplane. The effects of 
fluid viscosity, time lag, asymmetric forebody vortices, asymmetric 
leading-edge vortices, viscous flow separation, shock boundary layer 
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interaction, and other unsteady aerodynamic interferences should be 
somehow accounted for. At least, a theoretical method for predicting 
the important 0-derivatives for slender wings should be developed. 

It is shown that wing rock can be effectively suppressed if a roll 
damper SAS is incorporated to augment roll damping. However, more 
research should be done on high-a control law to provide a fast and 
reliable AFCS for future fighters. 
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Appendix A. Derivation of Approximate Analytic Solutions of the 
One-Degree-of-Freedom Model 

The Beecham-Titchener method (ref. 1) is basically an extension of 
the Krylov-Bogoliubov asymptotic method. This method is shown to be a 
combination of the averaging principle and the method of variation of 
parameters (ref. 42). For a damped linear system (ref. 1), the first- 
order approximation of the Beecham-Titchener method gives the exact 
solution. For such a system, the Krylov-Bogoliubov method yields an 
incorrect undamped natural frequency. Other examples (ref. 1) involving 
nonlinear terms are solved with excellent accuracy by using the Beecham- 
Titchener method. This method is also capable of obtaining the higher- 
order approximations. Therefore, this method is applied here. 

Equation (9) is rewritten as 

+ = L 0 + L s 6 a + sin a s h * + (L P o + sin “s + L p pl < t ) l ) + (A>1) 

A 

The solution of equation (A.l) is assumed to be 

4>(t) = a(t) cos v(t) (A. 2) 

a) and X are defined as 

0 ) = v (A. 3) 

X = a/a (A. 4) 
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where the superscript dot means a time derivative. In addition, let 

differentiation with respect to the amplitude a be denoted by a 

• • •• 

superscript prime. Then to, X, and a can be expressed as 

0 ) = a X (o' (A. 5) 


a X X' 


(A. 6) 


a = aX(aX'+X) 


(A. 7) 


On differentiating equation (A. 2) with respect to time, one obtains 

<(> = a(X cos v - a) sin v) (A. 8) 

<)> = a { [ X ( a X' + X) - <jJ^] cos v - X(auj' + 2u>) sin v} (A. 9) 

The Beecham-Titchener method assumes that uj and X do not vary 
greatly during one cycle of oscillation. In other words, a, X, w, and 
amplitude derivatives are interpreted as having fixed values and are 
assigned the values at midcycles. Equation (A.l) without the control 
deflection term is rewritten as 


4 » - L. 


+ sin o L„ a cos v + 
s 8 


[L 


'pO 


+ sin o L „ sgn (d>) a cos v 
s p8 


+ L sgn (<j>) a(Xcosv - oosinv)] a(Xcosv - msinv) 
PP 


= F (v) 


(A. 10) 
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Equations (A. 9) and (A. 10) are combined into 

F(v) = a[A(aX' + A) - oj^] cos v - aA(ao)’ + 2o>) sin v (A. 11) 

Regarding equation (A. 11) as a Fourier expansion for F(v), the Fourier 
coefficients can be determined to be 

o i 2 tt 

A(aA' + A) - a) = — / F(v) cos v dv (A. 12) 

** 0 

1 2tt 

A(au>' + 2oi) / F(v) sin v dv (A. 13) 

" a 0 

Equations (A. 12) and (A. 13) embody the averaging principle. They are 

similar in form to the Krylov-Bogoliubov first-order-approximation 

2 

equations, but include new terms in A and amplitude-derivatives. The 
amplitude derivatives are neglected in the first-order approximation. 
Mathematically, this is in accord with the method of variation of 
parameters (ref. 42). 

Before equations (A. 12) and (A. 13) can be evaluated, the following 
three groups of nonzero integrals have to be estimated in advance. 

Group 1 

2ir 

f sin v sin v dv = it 
0 

2ir 

/ cos v cos v dv = IT 
0 
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(A. 14) 


2 n 3 2 

f sin v sin v dv = yiT 

0 T 


2ir 


j cos J v cos v dv = 


-5-W 


2tt 


^ £ 2 I 

/ sin v cos v dv = ^ 

0 


7 1 


Group 2 


2ir 

f sgn (cj>) cos v dv = 4 
0 


2 ir 


2 

/ sgn ( 4 > ) sin v cos v dv 
0 


4 

3 


(A. 15) 


2 * 2 
/ sgn (<j>) cos ' 1 v cos v dv 

0 


8 

3 


Group 3 


2 ir 

/ sgn (<)>) sin v dv = -4 
0 


2it . 2 

/ sgn (<j>) sin 4 v sin v dv = 


(A. 16) 


2tt . 2 

/ sgn (<j>) cos v sin v dv = 
0 


4_ 

3 


After lengthy integration, equations (A. 12) and (A. 13) are reduced 


to 


A* - 


2 8 8 
a) “ sin a L„ + A L _ + aA sin a L „ + aA u L 

s 8 pO 3ir s pg 3 tt pp 


(A. 18) 
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(A. 19) 


2Xoj » a) L A + 4- a u sin a L . + 4— a(X^ + 2 oj^) L 

pO 3ir s p8 3tt pp 

All the aerodynamic derivatives in equations (A. 18) and (A. 19) should 
then be evaluated at some average conditions. 

Equation (A.l), a nonlinear second-order differential equation is 
then converted into two coupled nonlinear first-order differential 
equations, ft and A are defined as the circular frequency and amplitude 
for a steady-state wing rock. Limit-cycle characteristics are obtained 
by setting X to zero in equations (A. 18) and (A. 19). Equations (A. 18) 
and (A. 19) are then reduced to 


and 


■? 

-fi = sin a L„ 
s 8 


0 = ft L . + — A SI (sin a L „ + 2S2 L ) 
pO 3 tt s p8 pp 


(A. 20) 


(A. 21) 


The closed-form solutions of SI and A are 


S2 = (-sin ct L.) 

s 8 


1/2 


, qSb . „ ,1/2 

“ (“ i sin a s c , > 

l xx s *8 


(A. 22) 


A = - 


^ 0 _ 


4 sin a L . + 2S1 L 

s p6 PP 


3tt 

4 


'pO 


sin a C + rr— C 
s Z „ V Z 
pS PP 


(A. 23) 


It is noted that at limit-cycle conditions equations (A. 18) and (A. 19) 
have been reduced to the frequency and the amplitude equations 
respectively. 
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Appendix B. Derivation of Approximate Analytic Solutions of the 


Three-Degree-of-Freedom Model 

The equations of motion, equations (36), (37), (38), (40), and (41) 
are rewritten without control-deflection terms as 

8 - 5 0 + V + V * V * V (Ba) 

p • l 0 * V + <E p0 + Sb | 91 + £ PP |p|) p + V < B - 2 > 

' * V V * <\o + 5 pB 1 8 1 + "ppl p l> p + V <B - 3 > 

• • 


P = <$> (B.4) 

r = i|> (B.5) 


Solutions to the above equations are assumed to be 


0(t) = a(t) £(t) cos [ v(t) + e(t)J 


<j>(t) = a(t) cos v(t) 


t(t) = a(t) n(t) cos [v(t) + 6 ( t ) ] 


(B.6) 


(B.7) 


(B.8) 


Again, let m and X be defined as 

w = v (B.9) 

X = a/a 
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Also, all time derivatives are to be expressed in terras of amplitude 
derivatives as 

ui = a X ai ' 

• 

X = a X X* 

a = a X (a X’ + X) (B.10) 

• 

n = a n ' X 

n = a X ( aq "X + an ' X ' + n ' X ) 

• 

5 = ar x 

c = aX (ag "X + aC’X' + £'A) 

Because the oscillation frequencies of (j> , i|i, and B modes are assumed to 
be the same, e and 5 are virtually independent of time. To the first- 
order approximation, the amplitude derivatives are all neglected in 
equation (B.10). It follows that the time derivatives of equations 
(B.6) - (B.8) can be shown to be 


0 » a £ [(Xcose - oisine) cosv - (Xsine + oicose) sinv] (B.ll) 

• 

<j> » a(Xcosv - oisinv) (B.12) 

<j> ® a[(X 2 - ui 2 ) cos v - 2Xwsinv] (B.13) 

\p » an[(Xcos6 - uisin6)cosv - (Xsind + uicos<5)sinv] (B.14) 

2 2 

i|> => an{[X - m )cos 5 - 2Xoisin6]cosv 

- [(X 2 - m 2 )sinS + 2Xoi cosS]sinv} (B.15) 


Substituting equations (B.4) - (B.15) into (B.l) - (B.3) and averaging 
over one oscillation cycle according to Beecham-Titchener technique 
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yield the following first-order-approximation equations without control 
deflections as 


£(Xcose - cosine) a Y . £ cos e + Y X + Y n (Xcos6 - cosind) + Y, 

p p r $ 


£(Xsine + cocose) » Y. £ sin e + Y <o + Y n (Xsind + cocosd) 

p p r 


29 — — 8 — 8 — — 

X -co « L 0 £cose+L _X + ■=— L a a£X + ■=— L aXco+L r)(Xcosd-cosind) 
P pu pit pp jit pp r 


- - 4- 4-28-2 

2Xco»L„ £sine + L rt co + — L „ a£to + — L aX + ■=— L aco 
8 pO 3ir p8 3it pp 3n pp 


+ L n (Xsin6 + coc os6) 
r 


0 [ (X^-co^)cosd - 2Xcosind] » £ cose + N X +• 5 a£X 


+ 4— N aXco + N n(Xcosd - cosind) 
3ir pp r 


2 2 - 4 - 

rt[(X -co )sin<5 + 2Xcocos<5] » £sine + N _<o + r— N „ a£co 

P pu JTT pp 


+ 4— N aX^ + -f— N ato^ + N n (XsincS + cocosd) 
3tt pp 3 tt pp r 


(B.16) 


(B.17) 


(B.18) 


(B.19) 


(B.20) 


(B.21) 


During this integration process, a, £, n , X, co, e, and 6 assume fixed- 
values at raldcycles for slowly varying motions. In the meantime, | @ ( is 
approximated by |£<j)J because <j> and 8 modes are in phase in most cases. 
Again, all the aerodynamic derivatives should then be calculated at some 
average conditions. 
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Define G and H as the amplitude ratios of the 0 and 9 modes to 
the 4> mode at limit-cycle conditions. The phase angles, z and 6, are 
numerically negligible at limit-cycle conditions. In addition, X is 
zero for a steady-state wing rock. Equations (B.16) - (B.21) are thus 
reduced to 


0 = Y 0 G + Y 

p 9 

G = Y + Y H 

P r 

fi 2 = - L.G 

P 

0 = L n + 4~ L AG + 5— L AS1 + LH 
pO 3rr p0 3ir pp r 

-Hft 2 = -N G 

P 

0 = N„+4-N o AG+|-N Afl + NH 
pO 3-rr p0 3 tt pp r 


(B.22) 

(B.23) 

(B.24) 

(B.25) 

(B.26) 

(B.27) 


It is implied that the frequencies of 9, 0, and 9 modes are the 
same in equations (B.16) - (B.21). Because equations (B.22) and (B.26) 
are related to the frequencies of 0 and 9 modes respectively, they will 
be ignored in determining the effects of 0 - 9 amplitudes on 
the 9~aroplitude for a steady-state wing rock. Thus, equations (B.23), 
(B.24), (B.25), and (B.27) are reduced to 


G = Y + Y H 
P r 

= sin a - H cos 2 a + C + H C cos a 

3 3 2raV 2 Y p 2 mV 2 Y r 


(B.28) 


sr = - LG 

8 


= - G qSb 


^ \ ^ °n B 

I I * I 2 

XX zz xz 


(B.29) 
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A = - 


3n 


L n + HL 
pO r 


G L . + 2Q L 
PS PP 


3n 

4 


§ c £ + f c n 

po — -P-O- 


G(S C * « + f C n > + T- (g c £ + f c n > 

pg p6 pp pp 


where 


h = 


I C + I C 

zz £ xz n 

r r 

I C + I C 

xz l xx n 

r r 


g = 1 - h 


xz 


zz 


I I 

f = y!* - h) 

ZZ XX 


and 


H = - — [N . + 4“ A(G N + 2J2 N )] 
N P° 3u P 8 pp 

r 


I c. +1 c -4 -a[g(i C. +1 C )+^<I C. +1 c )] 

xz i xx n ^ 3 tt xz £ „ xx n „ V xz i xxn 

pO pO pB p8 pp pp 

(I C. + I C ) rne rv 
v xz Z xx n ' cos a 


(B.30) 


(B.31) 
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Appendix C. An Iterative Scheme for Calculating Average Aerodynamic 


Derivatives and Limit-Cycle Characteristics 

For a given angle of attack, it is necessary to calculate three 

values of C at different g and p in order to determine C and 
P pB 

by the VORSTAB code. The three values of are 
PP P 

called C, , C and r respectively. Due to the nature of 
pO pi p2 

nonlinear interaction between aerodynamics and dynamics, the most 

suitable g and p for calculating C , C and C are obtained 

^pO X pl p2 

through iterations. 

(J) and are defined as 


<f> = a cos v 

and 

• • 

<f> = a cos v - aui sin v = a(Acosv - cusinv) 


(C.l) 


(C.2) 


By applying the limit-cycle condition (X = 0) or for slow-varying 
amplitude with 4 => 0, equations (C.l) and (C.2) can be combined into 


? 2 
+ 2 = A 

52 


(C.3) 


or 


|i| = 52(A 2 - $V /2 


(C.4) 


where A and 52 are the limit-cycle amplitude of bank angle and circular 
frequency respectively. 
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Equation (C.3) represents the larger ellipse in the phase plane as 
shown in figure C.l. The smaller ellipse in figure C.l having one half 
of the limit-cycle amplitude of the larger ellipse can be defined by 


<J > 2 + 


*1 


1 A 2 

T A 


(C,5) 


From equation (C.4), the maximum roll rate is 


p 

max 



QA 


(C.6) 


The average bank angle $ ave and the average roll rate p ave are 
calculated from equation (C.5) to be 


ave 


i A 


and 


P « <fc 
ave ave 


2 “ ’ 2 <W 


(C.7) 


(C.8) 


The maximum and average sideslip angles, 3 „ and 8 , are obtained as 

max ave 


and 


max 


A sin a 

s 


0 = <b sin a = -i- A sin a 

ave ave s 2 s 


(C.9) 


(C.10) 


The dimensionless maximum and average roll rates are defined by 


and 


max 


°max 2V 


(C.ll) 
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^ave 13 ave 2V 


(C.12) 


Now C„ 


'po 


is calculated at g = Q and p = p 


. C„ , is calculated 
ave Z , 

Pi 


at 6 = 3 and p = 0. In actual calculation, this is taken to be 
ave 

at p = 0.05. C. is calculated at g = 8 and p = p 
K Z . ave K 'max 

p2 

Finally, C and C are calculated by 

A X 

pe pp 


c = _Rl _£0 

%8 ^ave 


(C.13) 


and 


C £ ~ 

_£2 l£l 


(C.14) 


PP 


max 


After C. , C» , and C. are obtained, the roll amplitude of wing 
>0 p8 pp 

rock based on the one-DOF model has been shown (eq. (18)) to be 


A = - 2l 

A f 4 


'p0 


, „ L bn _ 

sin a C. + 77- C„ 

s Z a V Z 

PS PP 


(C.15) 


Define 8^ and p^ as 


Bf = A f sin a g 


(C.16) 


and 


p f = 51A f b/2V 


( C . 1 7 ) 
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To start the iteration (figure C.2), a proper value of g is 

r r ave 

chosen, say 10 degrees. The corresponding p and p are calculated 

ave Mnax 

from equations (C.6) - (C.12). $1 is obtained from equation (16). 

C , C , and C are then properly calculated by using the VORSTAB 
p0 p8 PP 

code. The predicted C. , C p , C. , and C are used to 

g po P g pp 

find A.., 0 c ,and p_. If 0. and S (or p. and p ) do not agree within a 
f f f f max K max 

required accuracy, say 0.5 degrees, a second g is selected to be the 

nicix 

average of the first B max and g^. Two or three iterations are normally 
needed to get a final converged solution. 


In case that p is less than 0.05, C will be calculated 
ave Z . 

p° 

at g = 0 and p =0.05 and C. will be calculated at g = g and 
*ave Z „ ave 

P 2 

p = 0.10. In this case, only g, and g are iterated to be within the 
r max f max 

required accuracy. 
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Figure C.l. A Typical Phase Plane 




CONVERGED SOLUTIONS 
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Figure C.2. A Sketch for an Iterative Scheme. 



Appendix D. Energy Balance of a Steady-State Wing Rock for the One- 


Degree-of-Freedom Model 

The total aerodynamic energy exchanged over a limit-cycle oscillation 
is given by 


• 

AE = qSb / C (t) 4> ( t ) dt 
C(t) 


(D.l) 


or 


AE = qSb / C (<J>) d<j> (D.2) 

CC4> ) 

where C(t) and C(<j>) are the contours obtained by plotting versus t 
and <(> respectively. For a steady-state wing rock, C(<j>) is a closed curve 
over one oscillation cycle. Since the contours for positive and negative 
bank angles are theoretically the same, equation (D.2) can be simplified 
to be 


AE = 2qSb / C ($) d<)> (D.3) 

ABCBD 

where the integration path ABCBD is shown in figure D.l. The arrows in 
figure D.l indicate the direction of increasing time. 

The total aerodynamic roll damping coefficient is defined as 


C 


i 


Pt 


c o + c |0l + c |p| 

po x p6 x 


(D.4) 


The total rolling moment coefficient C 


i 


can now be rewritten as 
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p 


(D.5) 


C = C 


„ +■ C„ g + 

£ o £ b 


pt 


<j> derived in appendix C is given by 


\i\ - 0 (A 2 - <|) 2 ) 1/2 


(D.6) 


When <|> is increasing, e.g. along the path ABC, 


P 


SUA 2 - ♦ 2 ) 1/2 


(D.7) 


When 4> is decreasing, e.g. along the path CBD, 


P 


-a (A 2 


*V /2 


(D.8) 


Thus equation (D.3) can be rewritten as 


AE 


2 ’ Sb / ABCBD (0 * + C - > d » 


ABC 


'CBD 


(D.9) 


or 


A 

AE = 2qSb / (C - C ) d<fr (D.10) 

0 ABC CBD 


where 


'ABC 


\ + C £ 6 + C £ 

0 B pt 


bfi 2 .2.1/2 

2V (A " ♦ } 


CD. 11) 


and 
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(D.12) 


'CBD 


C £ + C * 6 
0 0 


~ bn /. 2 , 2 s 1/2 

C * 2V (A + 1 

Pt 


so that 


AE = 2qSb 2 | J"~ C ( - ^) U c d<*> 

V 0 pt 


,2 , 2 S 1/2 


(D. 13) 


The physical mechanism of a steady-state wing rock shows that the 

dynamic system is undamped over small bank angles and damped over larger 

bank angles. This can be explained by the dynamic loops in figure D.l 

with the help of equation (D.13). The loop ABD is clockwise and the area 

within this loop is positive. It results in a positive AE and implies 

that C is positive, or undamping. Thus, a clockwise loop is 
Pt 

dynamically destabilizing because the aerodynamic energy is added to the 

dynamic system. On the other hand, a counterclockwise loop, say the loop 

BCB, is dynamically stabilizing due to a negative C so that the 

Pt 

aerodynamic energy is consumed by the dynamic system. 

Substituting equations (D.4) and (D. 6 ) into (D.13) yields 


AE=2qSb 2 | / {[C 


0 


+sina C„ $■ 
l n s l a 

pO P0 


bfl 

lv 


(A 2 -tJ> 2 ) 1 / 2 J ( A 2 -<j>2) l/2 }d( j ) 


(D. 14) 


PP 


4 

where 


A 

r 

0 ~p 0 


/ C (A 2 -4> 2 ) 1 / 2 d4> = - 5 C {<KA 2 -ej > 2 ) 1/2 + A 2 sin " 1 A} 

1 pO ( 


= j C A 2 
4 % n 

pO 


(D. 15) , 
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(D. 16) 


r A , O w A 2 ,2 x 1/2 , 1 . _ rr , A 2 ,2x3, 1/2, A 

J sina C 4> ( A -<j> ) dt}> = — — - sina C {[(A -<j> ) ] } 

_ S X» _ j S 36 ^ _ 

0 pB pB 0 


-=- sin a C„ 
3 s l 


pg 


r A “ c 
o 2V 1 


PP 


< A 2 <{> 2 ) dcf) = || 


C £ { a 2 4> - | $ 3 } 

PP o 


_ i bn _ 

'3 V L Z 

PP 


(D.17) 


Finally, the aerodynamic energy exchanged over an ideal limit cycle is 


AE = 2qSb 2 §{1 C £ A 2 + isina g C £ A 3 + } H C £ A 3 } 
pO pB pp 


2 - , 2 .2 SI r 3Tr , . bSl „ ,, 

= ~ qSb A — {- — C + A (sina C + — C )} 

3 4 V l „ s l a V l 

pO P B PP 


2 - , 2 2 SI ,3tt 3tt _ , 

3 q Sb A v ( 4 C A - 4 C } 

pO pO 


= 0 


(D. 18) 


The intersection of stabilizing and destabilizing areas, i.e. the 

point B in figure D.l represents a turning point of roll damping 

because C is zero at that particular point. The critical bank and 
Pt 

sideslip angles, <j> and 8 , are defined as the angles of zero total 
c c 

aerodynamic roll damping. To find <j> , C. is set to zero to obtain 

C 

pt 

C - C - c k c l + § c (A 2 - * 2 ) 1/2 (D. 19) 

pt pO pB pp 

= 0 
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or 


l*. 


c * 0 + § c * <a2 - * 2 > 1/2 

__E° E 2 


-sin a C. 
s Z 


(D.20) 


pB 


Squaring both sides of equation (D.20) and rearranging give 


r <J> + s 6 + d = 0 

T c Y c 


(D.21) 


where 


P& PP 


s = 2 sin a C„ C„ 

s Z _ Z . 
pO pS 


j — c^- /bl a ri \ 2 

d ■ c * ‘ % A S > 

pO PP 


(D.22) 


Solving the quadratic equation (D.21) and comparing its solutions with the 
numerical solutions obtained from equation (D.20), <|> can be shown to be 


!<g 

1 , .2 
= "27 (-s + (s 

- 4rd) l/2 ) 



C. 

IU 

s 

pO 

2r - sin 

V* R 

PB 

i* c i 

1 f ( 2 
" 2? ( - S ‘ <S 

- 4rd) ) 


if C £ > 0 
PP 


if C £ =0 (D.23) 


i f C £ <0 
PP 
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Stabilizing Area 
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Figure D.l. A Typical Histogram of the Total Aerodynamic Rolling Moment Coefficient versus 

Dank Angle for One Limit Cycle of Wing Rock (The Arrows Indicate the Direction 
of Increasing Time). 



Appendix E. Suppression of Wing Rock with a Roll Damper 

A conventional roll damper is usually designed to quicken the roll 
response and augment roll damping. It is also a simple control device for 
suppressing wing rock. An ideal roll damper consists of a rate gyro which 
senses the airplane roll rate p and then feeds it back to actuate the 
lateral control unit 6^ in proportion to p but such as to oppose p. The 
lateral control deflection angle is then 


A6 


A 


Kp 


(E.l) 


where K is the system gain (control surface deflection per unit roll rate) 

of a roll damper. The corresponding increases in C^ , C^ , and due 

P P n p 


to A <5 „ are expressed as 
A 


4C Y ■ V • 46 A ■ 4C Y P 

5 a P 


4C t ’ % ' 45 a * 4C * p 

S A P 


AC = C • A6. = AC p 
n n r A n 

5 a p 


where 


P = P 2V 


AC = K — 
AC Y K b/2V 
P 


AC . = K b72? 
p 


(E.2) 

(E.3) 

(E.4) 

(E.5) 

(E.6) 

(E.7) 
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(E.8) 


iC „ " K b/fv 

P 


By replacing 5 , C , C. , and C in the equations of motion, equations 

A p pO n P o 

(25) - (27), by A6 a> C y + AC y , C £ + AC £ , and C n + AC n , the 

P P pO p pO p 

modified solutions at limit-cycle conditions are obtained to be 


A = - — 


Jn 

4 


g(C t + 4C ) + f(C„ + AC„ ) 
pO P pO P 


G(gC c + fC 

pb 


n > + r- (gc * 

P6 PP 


+ fC ) 
n 

PP 


(E.9) 


_ T zz C £ + l xz C n 

Q = -GqSb 2 

*xx *zz *xz 


(E. 10) 


G = sina - H cos 2 a + ^ Sb -r (C + AC ) + H — - b — C co s a (E.ll) 
S 3 2nw Y p p 2mV Z r s 


H= - 


I (C„ +AC )+I (C +AC )-4-A[G(I C„ +1 C )-f~(l C 0 +1 C )] 
xz l n l xx n - n 3ir xz £ . xx n . V xz £ xx n 
pO p pO P p8 p8 PP pp 

(I C„ + I C ) cos a 
xz £ xx n 
r r 


(E. 12) 


where 


h = 


I C + I C 

zz £ xz n 

r r 

I C + I C 

xz £ xx n 

r r 


g = 1 - h 


xz 


zz 


f - (Ih. - h) 

ZZ XX 


(E. 13) 


(E. 14) 
(E. 15) 
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Equations (E.9) - (E.12) are then used to calculate the system gain 

needed to alleviate the limit-cycle amplitude. A delta wing of 76° 

sweep with a vertical tail (ref. 29) is chosen to demonstrate the 

concept. Without a roll damper, C„ and AC. and the limit-cycle 

PO P 

amplitude are predicted respectively to be 0.0071 and -0.0410 and 50° 
at a of 20°. In contrast to that, the limit cycle can be eliminated 


A C Z 

completely by installing an ideal roll damper with K (K = ^ — 2.) less 

than -0.15°/rad/sec. Note that C « 0, C = 0.0046 (per degree) 

Y r Z j. 


and C = 0.0008 (per degree) calculated from wind tunnel data in 
reference 29 have been used in equations (E.6) - (E.8). The value of K 
seems small because the wing span of the wind-tunnel model is only 1 .63 
ft and the test velocity is 50.2 ft/sec. It is shown that the effect 
of AC Is negligible and augmenting roll damping via AC. is the most 

I JC 

p p 


efficient. Therefore, wing rock can be suppressed by using a roll 
damper to artificially augment roll damping. 
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Wing-Rock Buildup of an 80°-Delta Wing at a = 27° (ref. 
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Figure 4. Asymmetric Leading-Edge Vortices for an 80°-Delta Wing (ref. 8). 








<*, <*eg 

Figure 7. C» and Z 9 versus a for an 80°-De1ta Wing (ref. 8). 
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Figure 10. Rolling-Moment Hysteresis. 
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Figure 12. Discontinuous Aerodynamics with a 
Lumped-Time-Lag Concept (ref. 37). 
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Figure 13. Numerical Simulation of Wing Rock (ref. 38). 
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Figure 15 .Damping- in-Roll Derivatives of a Delta Wing with 80-Degree 


Sweep (X^ g =0.4c) at 0.1. 
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Figure 16.0antping-in-Ron Derivatives of a Delta Wing with 
76-Degree Sweep (X|-g=0.40c) at 


c-d- 
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Theory (V0R5TAB) 


C j> (Body Axis) 
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(Stability Axis)- (X CG =0.39c) 


Data of Forced-Oscillation Tests (ref. 48) 

+ c £* sin ° ( Bod y Axis) 

P - 8 

k = 0.20*0.12 , = + 2.5° 


Figure 17. Damping -in- Roll Derivatives of a Delta Wing with 74-Degree 
Sweep at Mj'O.Z . 
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Theory (VORSTAB) 

— - C* (Body Axis) 

P 

— (Stability Axis) (X C( ,=Q.34c) 

p s 

Data of Free-to-Damp Oscillation Tests (ref. 49) 

C» +C». sin a (Body Axis) 

*p b 

k = 0.13, 4<f> = + 30° 

--O — C» (Stability Axis) (X rr =0.25c) 

P s 

k = 0.14 - 0.10, A<(> = + 30° 

Figure 18.Damping-in-Ron Derivatives of a Delta Wing with 60-Degree 
Sweep at 1^=0. 1. 
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Figure 19 . Longitudinal Aerodynamic Characteristics of a Delta Wing 
of AR=0.705 with Sharp Leading Edges at M w =0.1. X(,q=0.4c; 
Re=1.4xl0 5 . 
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One -DOF Theory for 
Limit-Cycle Conditions 
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Figure 24. Nonlinear Roll Damping Coefficients for a Delta Wing of 80-Degree Sweep at 


OnerDOF Theory 
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Figure 25. The Histogram of the Total Aerodynamic Rolling Moment Coefficient versus Bank Angle for One 
Limit Cycle of Wing Rock at Angle of Attack of 27 Degrees. 


One~DOF Theory 
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Figure 26. The Histogram of the Total Aerodynamic Rolling Moment Coefficient versus Bank Angle for One 




0 I I ! 

One -DOF Theory 



Figure 27.The Critical Bank and Sideslip Angles for a Delta Wing 
of 80-Degree Sweep during Steady-State Wing Rock. 
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Figure 28. Lift Coefficient (Data Obtained at Re=696,000) and 
Steady-State Amplitude of Wing Rock for a Delta Wing 
(AR=1.0) with a Vertical Tail at M m =0.1. X CG =0.30c. 
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Figure 29. Lift Coefficient (Data Obtained at Re=812,000) and 

Steady-State Amplitude of Wing Rock for a Cropped Delta 
Wing of 76-Degree Sweep (AR=0.33)with a Vertical Tail 

at M =0.1. X--=0.27c. 


BE 





(b) Cropped Delta Wing with a Vertical Tail. 
S t /S=0.134. 

Figure 30. Comparison of the Onset a with Data at M^O.l. 

Moment-Center Location is at Forward Gravity- 
Center Position. 
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